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THE  PHYSICS  AND  THE  CHEMISTRY  OF  ETWAC  CODE 
(AN  EARLY  TIME  WET  AIR  CHEMISTRY  CODE  FOR 
DISTURBED  AIR  CONDUCTIVITY  CALCULATIONS) 

1.  INTRODUCTION 

When  air  is  bombarded  by  a  pulsed  beam  of  electrons  or  charged 
particles,  a  conducting  ionized  path  is  generated  which  will  subse¬ 
quently  decay  and  deionize.  The  conductivity  of  the  plasma  is  deter¬ 
mined  by  the  electron  density  in  the  ionized  region.  However,  the 
electron  density  or  the  deionization  of  the  plasma  depends  on  the 
atmospheric  ions  generated  as  a  result  of  the  energy  deposited  in  air 
by  the  bombarding  charged  particles. 

The  air  ions,  their  forms  and  their  paths  in  the  deionization 
scheme  depend  on  the  air  density,  various  kinetic  temperatures  and  the 
density  of  the  minor  species  present  in  air.  Minor  species  and  im¬ 
purities  (H^O,  0^,  C02,  NO,  N02,  N2Q,  etc.)  are  generally  present  in 
the  atmosphere.  Their  role  in  the  deionization  of  air  must  be 
assessed. 

In  this  report  we  present  the  physics  and  the  chemistry  of  a 
code  developed  to  depict  the  role  of  the  water  vapor  in  the  deioni¬ 
zation  of  air.  This  code  is  developed  to  provide  calculations  on  the 
time  history  of  the  electron  density,  at  a  point  in  space,  during  the 
passage  of  an  electron  beam  in  wet  air.  The  code  is  named  ETWAC 
(Early  Time  Wet  Air  Chemistry)  to  emphasize  its  application  to  the 
early  time  phase  of  the  disturbed  air. 

The  report  describes  the  deionization  of  the  dry  air  in  Section 

2.  In  Section  3  the  processes  which  arise  upon  the  addition  of  water 
vapor  to  air  are  presented  to  depict  the  role  of  H20  in  the  deioniza¬ 
tion  processes.  In  Section  4  the  rate  coefficients  currently  used  in 
Manuscript  submitted  April  10,  1981. 
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the  code  are  presented  while  in  Section  5  a  discussion  is  given  for 
the  calculation  of  various  relevant  temperatures. 

II.  THE  CHEMISTRY  OF  DRY  AIR 

The  ionization  of  the  air  species  (N^  and  0^)  by  the  electron 
beam  arise  from  the  collisional  ionization,  ionization  due  to  the 
Bremsstrahlung  radiation  and  the  avalanche  ionization  caused  by  the 
intense  electric  field  associated  with  the  electron  beam.  Regardless 

■f  ^  *4* 

of  the  ionization  sources,  the  following  ions:  N  ,  0^,  and  0  ,  are 

generated  directly  in  the  dry  air  along  with  the  free  electrons. 

The  deionization  (i.e.  the  disappearance  of  the  charges,  positive 
and  negative)  proceeds  in  two  distinct  manners:  1)  by  the  electron- 
ion  recombinations  and  2)  by  the  ion-ion  recombinations.  These  pro¬ 
cesses^,  however,  are  complimented  by  various  atomic  processes  which 
result  in  the  generation  of  other  ions  not  produced  originally  in  air. 

The  electron-ion  recombinations  proceed  through  the  dissociative 
and  the  three-body  recombinations.  The  dissociative  recombinations 
are: 


e  +  N* 

-► 

N  +  N 

(1) 

e  +  N* 

-> 

N  +  N, 

(2) 

e  +  h’t 

4 

-*■ 

V  N2 

(3) 

e  +  0^ 

-*• 

0+0 

(M 

e  +  o'!" 

k 

-*■ 

°2  +  °2 

(5) 

e  +  N0+ 

N  +  0 

(6) 

2 


e  +  0+  H 


°2  +  N2 


(7) 


The  three-body  recombinations,  on  the  other  hand,  proceed  according  to 


e  +  e  +  A  ■+•  e  +  A 
n  n 


and 


e  +  M  +  A 


M  +  A 


n 


(8) 


(9) 


Where  the  third-body,  M,  in  equation  (9)  is  a  neutral  species  and  A 

n 

indicates  an  atomic  ion  (n=l),  a  molecular  ion  (n=2)  or  a  cluster  ion 
(n“3  and  4) .  Recombination  also  proceeds  through  the  radiative  re¬ 
combination,  which  can  be  ignored  in  regions  of  interest  to  our 
problem. 


ions 


..+ 

ions, 

»  N+, 

0*  and  which 

appear  in  equations 

i  (7)  are 

formed  by  the  following 

association  react- 

N+  +  N2 

+  N2 

CM 

+ 

V* 

+ 

(10) 

N2  +  N2 

+  N2 

-  N+  +  N2 

(ID 

°2  +  02 

+  02 

-  0t+°2 

(12) 

CM 

+ 

V 

+  N2 

-*■  02.N2  +  N2 

(13) 

hand,  is 

formed  by 

0+ 

+  N2 

N0+  +  N 

(14) 

N+ 

+  02 

"*■  N0+  +  0 

(15) 

and  by  various  charge  exchange  processes  of  positive  ions  with  NO. 
However,  for  early  time  phase  of  the  beam  air  interaction  one  may 
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ignore  NO  . 

The  ion-ion  recombination,  also  called  the  mutual  neutralization 

can  be  illustrated  by  the  following  two  reactions 

0~  +  A+  ■+  0„  +  A  (16) 

z  n  z  n 

O'  +  A+  h-  0  +  A  (17) 

n  n 

Mutual  neutralization  also  proceeds  through  the  three-body  neutrali¬ 
zation  reaction 

(O', Or)  +  A+  +  M  -*■  (0,0_)  +  A  +  M  (18) 

z  n  z  n 

where  the  negative  ions  0^  and  0  are  produced  by  the  three-body 

e  +  °2  +  °2  -  0*  +  °2  (19) 

e  +  o,  +  N2  -  0~  +  N2  (20) 

and  by  the  dissociative  attachments 


e  +  02  0  +  0 


respectively. 


The  positive  ions,  however,  undergo  various  charge  exchange  and 
ion  atom  interchange  processes  which  are 


N2  +  °2  02  +  N2 


N+  +  02  -*■  0*  +  N 


0  +  o2  -  o2  +  0 

N"!"  +  0o  -*■  ot  +  2N„ 
4  2  2  2 


N*  +  02  N0+  +  N, 


(27) 


A  schematic  diagram  shown  in  Fig.  1  illustrates  the  flow  of  the 
positive  charge  in  the  deionization  scheme  of  the  disturbed  dry  air. 
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III.  THE  CHEMISTRY  OF  WET  AIR 

The  presence  of  water  in  air  introduces  further  complications  in  the 

disturbed  air  chemistry.  Various  positive  clustered  ions  are  formed  which 

terminate  in  the  form  of  a  hydrated  hydronium,  H.0+. (H„0)  .  These  ions 

o  Z  n 

2 

have  been  observed  by  Narcisi  and  Bailey  in  the  D-region  of  the  ionosphere. 

The  reaction  paths  leading  towards  these  clusters  starts  with  0^  as  a  pre- 

,3,4 


cursor  ion  and  is  well  understood" 


and  have  been  utilized  in  the  mod¬ 


eling 


1,5 


of  the  disturbed  atmosphere  as  a  result  of  nuclear  detonation. 


The  hydrated  hydronium  generally  recombine  with  the  free  electrons, 
through  the  dissociative  recombination,  at  a  much  faster^  rate  compared  to 
the  lighter  air  ions. 

The  pertinent  reactions  which  lead  to  the  generation  of  various  hy¬ 
drated  ions  are  as  follows.  The  positive  ions,  N+,  N*  N*,  N*  and  0+  react 

Z  J  *4 

with  H^O  according  to  the  following  reactions 


A 

A 

A 


+  h2o 

-V 

H20+  +  N 

(29) 

+  H20 

“► 

h2o+  +  n2 

(30) 

+  H20 

N2H+  +  OH 

(31) 

+  H20 

-V 

h2no+  +  n2 

(32) 

+  H20 

H20+  +  2N2 

(33) 

+  H20 

-4 

h2o+  +  0 

(34) 

Reactions  (29)  through  (34)  obviously  produce  predominantly  H^O  which  in 
turn  undergoes  the  following  reactions 

H20+  +  °2  °2  +  H2° 


i 


(33) 


OH  +  H 


(36) 


H90+  +  e 

H  0+  +  H90  -  H30+  +  OH  (37) 

Reaction  (37)  generates  the  hydronium,  H^o"*", which  in  turn  generates 
once  hydrated  hydronium  according  to 

H30+  +  H20  +  M  -  H30+-H20  +  M  (38) 

The  following  reactions,  an  association  and  two  switching,  result  in 
the  production  of  0*-H20 


09  +  H,0  +  M 

-*•  02'H2° 

+ 

M 

(39) 

°t  +  h2o 

-*■  o2.h2o 

+ 

M 

(40) 

02-N2  +  H20 

"  °2-H2° 

+ 

N2 

(41) 

"f  + 

On  the  other  hand,  02*H90  reacts  with  H20  to  produce  H^O  -OH  which 

in  turn  produces  H30+-H90  through  a  switching  reaction  with  H20.  This 

indicates  two  routes  for  the  production  of  H30  *^0,  one  initiated  by 

* 4 ~  *4"  "4* 

H30  ,  and  the  other  by  02<H20  where  0 2  and  0^  are  the  precursor  ions. 

Other  hydrates  are  generated  according  to 

H30+- (H90)n  +  H20  +  M  h.  H30+- (H20)n+1+  M  (42) 

All  these  new  clusters  recombine  with  the  free  electrons  through 
the  dissociative  recombination  process  to  produce  neutral  products. 

A  schematic  diagram  shown  in  Fig.  2  illustrates  the  flow  of  the 
positive  charge  in  a  disturbed  wet  air.  This  diagram  clearly  shows 
more  complication  compared  to  dry  air  scheme  presented  in  Figure  1. 
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+ 
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Fig.  2  —  Disturbed  wet  air  positive  ion  flow  charts 
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IV.  REACTION  RATES 


The  relevant  reaction  rates  pertinent  to  the  early  time  wet  air  chem¬ 
istry  code  is  given  as  a  function  of  temperature,  if  the  temperature  depen¬ 
dence  is  known;  otherwise  the  room  temperature  values  are  presented.  A 
review  of  these  and  other  rates  is  to  be  reported  elsewhere7. 

Table  1.  LIST  OF  REACTIONS  AND  THEIR  COEFFICIENTS 


Reaction 

Rate  Coefficient 

Reference 

N+  +  e  N  +  N 

,  -j  ,H8  rv  T0-39 

4. 3  x  10  (T  ) 
e 

8 

ot  +  e  -►  0  +  0 

1.5  x  108  (T  I0,7, 

T  <0.1 

8 

2 

e 

e  - 

2.1  x  108  CT  )0,5, 

T  >0.1 

8 

N*  +  e  -  N2  +  N2 

e 

3.4  x  108  CT  )1-1 
e 

e 

9,  10 

N ^  +  e  -*•  N  +  N2 

1.75  x  108CT  )1,0 
e 

6,10 

°4  +  e  °2  +  °2 

3.4  x  108  (T  T1*0 
e 

6,10 

H30+  +  e  -*-H  +  H20 

3.2  x  107  (Tj1*0  , 

Te  <  0.86 

6,10 

2.5  x  108  (T  )1-2  , 

0.26  <  T 

^2.15  11 

e 

e 

3.05  x  108(T  J1,43 

2.15  <  T 

<8.6  11 

e 

-  e 

H20+  +  e+  OH  +  H 

2.7  x  108  (T  )0’5 
e 

Estimated 

H30+’H2°  +  e 

1.6  x  106  (T  y0,15 
e 

6,10 

2H20  +  H 

H30+.  (H20)2  +  e  -> 

4.2  x  106  (Te)0,05 

6,10 

3H?0  +  H 
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Table  1.  (.continued) 


O*  .H20  +  e  -*■ 

7.2  x  107  (T  )0,2 
e 

6 

o2  +  h2o 

H30+-0H  +  e 

9.6  x  10?(T  T0,2 
e 

6 

h2o  +  H  +  OH 

e  +  M  +  (A+,A2)  ■* 

-31  -2  5 

5.9  x  10  (T  )  D 

e 

6 

(A,A2)  +  M 

e  +  e  +  (A+,A2) 

4.3  x  1027  (T  )4-5 
e 

6 

e  +  (A,A2) 

N2+2N2 

Nt  +  N2 

-35  -3  9 

3.1  x  10J^  (T  )‘i*y 
a 

1.3  x  1030  (T  )1,0 

cl 

12,13 

10 

N+  +  2N2  -f  N*  +  N2 

4.5  x  1031  (T  )1,0 
a 

10 

°2  +  202  ^  °4  +  °2 

1.9  x  1035  (T  )3,2 
a 

14 

02  +  M  +  H20  -»■ 

.  ,t31,^  r-2.0 

1.75  x  10  (I  ) 
a 

m=n2 

10 

02-H20  +  M 

1.43  x  1031(T  )2,0 
a 

m=o2 

10 

h3o+  +  h2o  +  M  -*■ 

2.1  x  1030  (T  )2,0 
a 

m=n2 

10 

H30+-H20  +  M 

1.4  x  1033  (Ta)4,0 

M=02 

10 

H30+-H20  +  H20  +  m  -*■ 

1.4  x  1030  (T  )2,0 
a 

m=m2 

10 

H30+. (H20)2  +  M 

7.8  x  1033  (T  )4,0 
a 

M=02 

10 

N+  +  02-  02  +  N 

2.8  x  1010 

T  <  0.39 

a  - 

10,15 
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Table  1.  continued 


5.3 

x  10“10 

(T  9  *  5  7 

T  >  0.39 

a 

a 

N2  + 

°2 

-  o2++n2 

2.7 

x  1012 

(T  )°-8 
a 

,  T  <  0.3 

a  — 

15 

4.2 

x  10U 

(T 

,  T  >0.3 

15 

a 

a 

Nt  + 

0o 

-*■  N0+  +  N  +0 

5.1 

x  10U 

16 

3 

2 

2 

and  N02  +  ^2 

8.8 

x  1012 

(I)0'52 

a 

17 

Nt  + 

°2 

-  0*  +  2N2 

2.5 

x  1010 

16 

0+  + 

Oo 

tm 

-*•  02  +  0 

4.6 

,n12 
x  10 

(T  ' )°‘4 

a 

t  <  .155 
a  — 

18 

1.0 

x  1010 

(T  P 

a 

T  >  .155 
a 

N+  + 

h2o 

H20+  +  N 

2.8 

x  10"9 

16 

N4  + 

H20 

H20+  +  2N2 

3.0 

x  109 

16 

0+  + 

H2° 

-  h2o+  +  n2 

3.2 

x  109 

16 

N2  + 

h2° 

-*■  h2o+  +  N2 

2.8 

-9 

x  10 

16 

and 

n!h  +  OH 

a=0 

.82 

Z 

b= 

.18 

h2o 

-  H2N0+  +  N2 

3.3 

■  xlO 
x  10 

16 

0+  + 

N2 

-*■  N0+  +  N 

7.5 

x  1010 

T  <0.1 
a  ~ 

15, 

3.2 

x  1011 

(T  P8 

T  >0.1 

a 

a 

T 

=  T 

a 

V 

°I  + 

H2° 

*  °2'H2° 

1.5 

-9 

x  10 

19 

12 


Table  1.  Continued 


°2‘H2°  +  H20  -  H30+-0H  + 

°2  1-2  x  109 

19 

^  HO  +  OH  + 

J  2 

,  a=0. 83 

b*0. 17 

H  0+-0H  +  H  0  -»• 

3  2 

1.4  x  109 

19 

H30  -H90  +  OH 

h20+  +  H2°  **  H30+  +  OH 

1.8  x  109 

10 

H2°+  +  °2  ^  °2  +  H2° 

2.0  x  1010 

10 

O+O  -»■  0  +  0 

4.3  x  108  [T  I0,5 

6,10 

4*  — 

N  +  0  -  N  +  0 

4.1  x  108  [TaJ0*5 

6,10 

°2  +  °~  ^  02  +  0 

1.6  x  108  [T  I0,5 
a 

6,10 

°2  +  °2  "  °2  +  °2 

6.6  x  108  [T  I0,5 
a. 

6,10 

N2  +  °2  "  N2  +  °2 

2.5  x  108  [T  T0,5 
a 

6,10 

X+  +  Y~  -*■  X  +  Y 

1.6  x  108  [T  T0'5 
a 

20 

4  — 

A+B+M-v 

A  +  B  +  M 

2.9  x  1029(T  J2,5 
a 

6 

+  — 

X  +  y  +  M  -*■ 

1.0  X  1029[T  T2,5 

1  aJ 

6 

X  +  Y  +  M 

e  +  o2  +  02+yo2 

3.6  x  1031(^-)  Exp (-^£12) 
c  e 

21 

e  +  02  +  N2  -*■  0~  +  N2 

1.0  x  1031 

21 

e  +  02  +  H20  -*»  02  +  H20 

1.4  x  1029 

21 
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Table  1.  Continued 


e  +  CL  -  0~  +  0  1.5xlO"i0(T  J1*0  Exp  (-5/T  )  22 

2  e  e 

T  and  T  are  in  eV 
e  a 

V.  SPECIES  AND  TEMPERATURES 

The  early  time  disturbed  wet  air  chemistry  code  solves  a  set  of  rate 

equations  to  calculate  the  time  histories  of:  N+,  N^,  N+,  0+,  0^,  0^, 

Ho0+,  H,0+-0H,  H  0+*Ho0,  0o  and  0.  On  the  other  hand,  the  electron 

3  3  22  32  2 

density  is  calculated  via  the  charge  neutrality  condition.  In  addition  to 
the  above  species,  the  electron  temperature,  T^,  the  N^  vibrational  temp¬ 
erature,  Tv, and  the  neutral  species  temperature,  T  ,  are  calculated.  The  set 
of  the  rate  equation  for  species  and  temperatures  are  solved  by  a  computer 
program  devised  by  D.  Strickland  of  Beers  Associates. 

The  calculation  of  T  ,  T  ,  and  T  are  necessary  because  (see  Table  1, 

e  v  a 

Section  4)  various  relevant  deionization  reactions  are  temperature  depen¬ 
dent. 

The  calculations  of  these  temperatures  requires  data  on  all  appropriate 
23 

cross  sections  for  elastic  and  inelastic  processes  in  order  to  account  for 
the  energy  loss  by  electrons  in  air.  These  cross  sections  are  then  utilized 
along  with  the  electron  velocity  distribution  to  obtain  the  relevant  rate 

coefficients.  The  current  rate  coefficients  used  in  this  code  are  based 
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on  a  previous  set  of  cross  sections  ,  some  of  which  have  been  modified 

to  the  most  current  data.  The  calculation  of  the  vibrational  temperature 

and  the  energy  stored  in  the  ^  vibrational  mode  follows  a  method  reported 
2  6 

previously.  The  neutral  temperature,  on  the  other  hand,  is  calculated 
ab  initio,  where  the  heating  sources  are  the  elastic  electron  neutral,  elec¬ 
tron  ion  collisions,  dissociative  recombination,  charge  excharge  and  the 
arrangement  processes. 


The  ohmic  heating  of  the  thermal  electrons  is  considered  in  the 


code  by  solving  a  simple  circuit  equation.  This  circuit  equation  re¬ 
lates  the  net  current  I  ,  to  the  induced  electric  field,  E  ,  via, 

n  z 


L 


dl 
_ _ n 

dt 


E 

z 


Where  L  is  the  inductance  per  unit  length  and  is  defined  as 


Where  r,  is  the  beam  current  radius  and  R  >  r.  is  the  radius  where  con- 
b  b 

Q 

ductivity  in  the  plasma  sheath  falls  to  -  .  The  net  current,  on  the 


other  hand,  is 


I  =  I.  +  nr,  C  E 
n  b  b  z 


a  is  the  conductivity  and  I  is  the  beam  current. 
The  conductivity  c  is  defined  as 

c 


2  » 
e  Ne 


m  Vm 


Where  is  the  electron  momentum  transfer  collision  frequency  and  is 
the  Siam  of  the  electron  neutral  and  electron  ion  collision  frequencies. 
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